Tsallis mapping in growing complex networks with fitness 
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Abstract 

We introduce Tsallis mapping in Bianconi-Barabasi (B-B) fitness model of growing networks. 
This mapping addresses the dynamical behavior of the fitness model within the framework of nonex- 
tensive statistics mechanics, which is characterized by a dimensionless nonextensivity parameter 
q. It is found that this new phenomenological parameter plays an important role in the evolution 
of networks: the underlying evolving networks may undergo a different phases depending on the q 
exponents, comparing to the original B-B fitness model, and the corresponding critical transition 
temperature Tq is also identified. 
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I. INTRODUCTION 



The seminal work by Watts and Strogatz on small- world networks and that by 
Barabasi and Albert on scale-free networks (i.e., B-A networks) |2|, have pioneered a wave 
of research activities on the complex networks in last decade (see, for instance, recent review 
articles {3 — 5 1 , books js, 7|, and references therein). Research toward understanding the new 
critical phenomena in complex networks, e.g., the structural phase transitions in networks, 
is also a rapidly developing area js]. It is found that the condensation phase transitions 



occur m many comp 
(B-B) fitness model 



ex network models U|22j. Particularly, in Ref. 



□ 



Bianconi-Barabasi 



12] was mapped into an equilibrium Bose gas, a Bose-Einstein conden- 



sation (BEC) behavior was obtained, and the critical point for condensation was identified 
by mean-field arguments. There the well-known B-A networks are generalized by intro- 
ducing the fitness parameter, rji. At each time step, a new node with m links is added to the 
nodes that already presented in the network. The fitness r)i, accounting for the competition 
for links in the network, is chosen randomly from a distribution. The probability H that a 
new node connects one of its m links to a node i depends on the number of links, fc;(t), and 
on its fitness 77,, is 

n, = (i) 



In order to map the fitness model into Bose gases, one assigns an "energy" e^, introduced 
randomly from a distribution g(e), for the zth node, and the fitness of the node, r]i, is defined 
as 



(2) 



where the parameter (3, can be identified as inverse "temperature" (3 = 1/T. This grow- 
ing network can be solved in a mean-field approximation and a "chemical potential" fi is 
determined by the self-consistent equation 

I(P,H) = J deg(e)n(e) = l, (3) 

where the Bose occupation number n(e) = 1/ (e^^ e ~^ — 1) is the number of links attached by 
the preferential attachment mechanism to nodes with "energy" e. A BEC phase transition 
occurs in the networks if I ((3,0) < 1, i.e., the self-consistent equation (JSJ) has no solution 
(in such case the self-consistent approach fails). A necessary condition for condensation is 



2 



that the distribution g{e) — > for e — > 0. For this type distributions there exists a critical 
temperature T c = 1/(3 C such that I ((3, 0) < 1 for T < T c . 

However, it should be noted that BE statistics appearing in B-B growing network is only 
formally equivalent to that in ideal Bose gases. This is due to the fact that the evolution of 
B-B network is in essential irreversible and out-of-equilibrium, while the BEC of ideal Bose 
gases is a temperature-driven equilibrium process. On the other hand, it was argued [23] 
a wide variety of complex networks have self-similar, fractal property. Under such circum- 
stances, a particularly useful and appropriate framework to describe the d yna mical evolution 



of the network is so-called Tsallis nonextensive statistical mechanics 



2J-|26|. It is a certain 



generalization of the usual Boltzmann-Gibbs (BG) formulation. For the applications of the 



27|, 



28| and 



formulation in many fields, including complex networks, one may see, e.g., Refs. 
references therein. In present work, we consider the possibility of embedding the B-B fitness 
model in the framework of Tsallis formulation of the non-extensive statistical mechanics. 
The paper is organized as follows: in Sec. HI1 we briefly review Tsallis formulation of the 
nonextensive statistical mechanics, and introduce Tsallis mapping in the B-B fitness model 
of the complex network. We show that the corresponding rate equation has a power-law 
solution under some approximation within the high temperature regime, in thermodynamic 
limit. For general situations without the thermodynamic limit, the numerical simulations 
are need. The main simulation results are shown and discussed in Sec. IIHI Finally, we 
summarize and draw our conclusions. 



II. TSALLIS MAPPING IN B-B FITNESS NETWORKS 

In Tsallis formulation of non-extensive statistical mechanics, the central equation is so- 
called the nonextensive entropy, which is defined as 

^=A(!-E^ ( 4 ) 

i 

where Pi is the probabilities for the ensemble to be in the state i. The non-extensive 
parameter q describes the degree of deviations of a physical quantity from extensive and 
therefore may be regarded as a new phenomenological parameter of the model used, different 
values of q — 1 quantifies the departure from the BG limit. The standard results of BG 
statistical mechanics are recovered in the limit q — > 1. The application of the non-extensive 
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statistical mechanics to Boson gases predicts a g-deformed BE statistics 251 ]: 



= ^ 7 , (5) 

-9 



where is the occupation number of an energy level £ and the chemical potential. The 
parameter B i s the inverse temperature. A g-exponential function is introduced in Eqn. (J5]) 



as follows 



29|: 



e « = [! + (! -^/M ( 6 ) 



if the argument 1 + (1 — q)z is positive, e z q equals to zero otherwise. Its inverse function, the 
g-logarithm, for any positive real number z G M + , is 

z 1 ~ q — 1 

ln g ^^--. (7) 

Both g-logarithm and g-exponential functions go back to the usual logarithm and exponential 
functions, respectively, when one takes q — > 1 limit. 

In order to map B-B fitness model to a boson gas with Tsallis statistics, under the self- 
consistency condition, we assign the ith node with a fitness, r^, by 

Vi = e" fe = (l + (1 - g)(-^)) V(1_<?) , (8) 

and its inverse is 

1 1 ri l ~ q - 1 

£i = -^ln,r fi = --^— , (9) 

where follows the energy level distribution g(e) = Ce e and C is a normalization factor. 
A necessary condition for condensation is that when e — > 0, g(e) — > consistently. The 
resulting rate equation which describes the degree of connectivity at time t; of the ith node 
(energy level), ki(Ei,t,ti), is then: 

dt tZ q t ' 1 j 



where 

Z3 = 

mt 



is the g-deformed partition function. Following the standard statistical mechanics 30], the 
definition of fugacity z is given, in a consistent way, by g-deformed form, z~ l = e~^\ in 
which fi is the chemical potential. 
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The phase structure of this B-B fitness model with Tsallis mapping is now controlled by 
two parameters, i.e., the inverse temperature /3 and the nonextensive parameter q, instead of 
(only) one parameter B in the original B-B fitness model. These two parameters both affect 
the phase transition of the network. A few limiting situations are worthy of being discussed 
before further detailed numerical simulations. First, in q — > 1 limit, the B-B model with 
Tsallis mapping goes back to the original one with Bose gas mapping, due to the fact that 
Tsallis g-deformed expressions approach to the usual ones, and all the known results of the B- 
B model are recovered as they should. Henceforth, the scale-free phase in original B-B fitness 



model is trivially recovered when all nodes have the same fitness, i.e., g(e) ~ 5(e) [2|, 111 ]- 
Second, with q — > 1 limit and the thermodynamic limit t — > oo as well, we find that the rate 
equation ffTU]) has a power law solution under some appropriate approximation, within high 
temperature region, in which B is a small quantity. This ph matter of fact, is similar 

to the ht-get-rich (FGR) phase in the orlglnal B-B fitn » m0 del Q. 
To see the latter, we assume that the solution takes the following form: 

kiieM)^ - , (12) 



A. 

where f q (s) is some unknown power- law exponent and to be determined later. For simplicity 
we will write f q (s) = f q from now on. Taking the energy distribution of the nodes g(e) into 
account, the partition function becomes 

Zf ~ J deg{e) j' ^V B A*(e, t, t ), (13) 

here we use Zl stand for the partition function averaging over g(e). Now inserting hi (in 
Eqn.([12l)) into Zf, one obtains 



Zl - I deg{e) / 

't 1 -^ - 1 



deg(e)e-P £ t fq 



-i 



= J deg(e)e^~\ln fq t) 

~ z- x {l-t-V-M), (14) 
where ln/ 8 1 is the /^-logarithm as introduced in Eqn. (J7J) and the inverse fugacity z^ 1 is 



e, 



-/3e 



z- 1 = I deg{e)^j. (15) 
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Note that in Eqn. (|14p . 1 — f q > 0, such that in thermodynamic limit t — > oo, t"^ 1 "^' — > 
0, and hence Zl ~ z -1 in the corresponding limit. This makes one able to determine 
aforementioned power-law exponent f q . Note that in high temperature region, (3 is small, 
and that q is very close to unity (such that q — 1 is small as well), so we can take the 
approximation e~ /3e /e~ /3/i ~ e q ^ [31], as a result, the exponent f q of the power-law 
solution in Eqn. ( fl2l) is, 

f q = e;n<-»\ (16) 



in q — > 1 limit, this exponent goes back to the one in the original B-B fitness model 
The fugacity becomes 



z- 1 = lim Z* = I deg{e) % (17) 

J 1 - Jq\£) 

and the chemical potential \x is the solution of the following equation, 

I q {fi,li) = I deg(e) ^ = 1, (18) 

which indicates a mapping to Bose gas with Tsallis statistics, and the corresponding oc- 
cupation number of a energy level e is n(e) = (e^ e ^ — l) -1 . This just follows Tsallis 
statistics. 

For general situations without taking into account above limits, one has to resort to 
numerical simulations for current model with Tsallis mapping. 



III. NUMERICAL SIMULATION RESULTS AND DISCUSSIONS 

In the original B-B fitness model, the signal of the condensation phase transition lies 
on the change of sign of the chemical potential /i, i.e., when \x experiences a change from 
negative value to positive value within some regime, that change identifies the critical point 
of the corresponding phase transition. Similarly, we use the same criterion for the phase 
transition in the B-B fitness network with Tsallis mapping. This is done by numerically 
computation of the chemical potential /i, without the thermodynamic limit, according to 
the following way: 

where ln g Zf is the g-logarithm of Z^ which is defined in Eqn. (II ip . Note that in above 
numerical simulations, the computation of fi strictly follows its definition and has no ap- 
proximations in them. 



The "phase diagrams" of the B-B fitness network with Tsallis mapping are demonstrated 
by our numerical simulations and are shown in Fig. [1] - [3j In both figures, the size of the 
evolving network, or the time steps, is taken to be t = 10 3 and we fix m = 2. In Fig. [TJ 
the absolute value of the chemical potential /i, which is defined in Eqn. ( TTTTj) . is numerically 
computed with an average over 300 runs, using the energy level distribution g(e) = 2e for 
energies e G (0, 1), i.e., 9 = 1 and C = 2, and for q = 1.0, 1.1, 1.3, and 1.5, respectively. 
As expected, for q = 1.0, the phases emerge in the system are coincident with those in the 
original B-B fitness model: the FGR phase and the BEC phase (the solid line in Fig. [1]), 
and the same transition temperature is located at Tbe ~ 0.8 (the vertical dash-line on the 
right). However, when the value of parameter q deviates from 1, e.g., for q = 1.1, 1.3, 
and/or 1.5, the phases correspondingly change depending on both the specific value of q 
and the temperature they take (see the dot-dashed line in Fig. [[]). Also, it can be seen from 
the figure that the transition temperature T c 9 decreases and T c 9 ~ 0.67 (the vertical dashed 
line on the left in Fig. [T|). In addition, in T < T c 9 regime, the whole "Tsallis condensate" 
curve deviates from the BEC curve. On the other hand, it is no surprise that the effect 
of q parameter disappears gradually in the high temperature limit. The effect of different 
energy level distribution is shown in Fig. [21 with 9 = 1/2. We also plot the fraction of the 
total number of links connected to the most connected (or "winner") node, k max /mt, as a 
function of T (lower panel in Fig. [T]and Fig. [2]). 

We also plot in Fig. J3] the connectivity degree distribution of the B-B fitness network 
with Tsallis mapping, i.e., the probability distribution P{k) vs. the connectivity degree k, 
for different values of parameter q and temperatures T. In the figure, P(k) is plotted for 
g(e) ~ e, with nonextensive parameter q = 1.0, 1.3, and 1.5, and with temperature T = 0.08, 
T = 0.4, T = 0.8, and T = 4.0, respectively. Since the critical temperature is Tc = 0.67, 
in the latter cases (i.e., T = 0.8 and T = 4.0), the temperature T > Tc- As we expected, 
the Tsallis mapping gives the power-law degree distribution for T > Tc, P{k) ~ A; -7 , with 
7 ~ f q (see the lower panels in FigJSJ in which T = 0.8 and T = 4.0, respectively). Note 
that, first, the nonextensive parameter q controls the slope of P(k) at given temperature; 
second, data points with q = 1 (open squares with blue color in Fig. |3]) in fact correspond 
to the FGR phase in the original B-B model for T > Tc- When the temperature decreases, 
for example, at the critical point Tc ~ 0.67, the chemical potential changes its sign from 
positive to negative, and this correspondes to the situation that a condensate phase starts to 
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FIG. 1: Numerical simulations of the fitness model with Tsallis mapping for growing complex 
networks. An energy level distribution g(e) ~ e is taken in the figure. We set evolving times (system 
size) to be t = 10 3 , and averaged over 300 runs. Upper Panel: Numerically calculated chemical 
potential using the energies distribution with e £ (0, 1). The temperatures at which fi changes sign 
correspond to the sharp drop in \fj,\ on the figure, and identify the different critical temperatures 
Tc (two vertical dashed lines in the figure) for Tsallis condensation phase with different g-value: 
q = 1.0 (solid line), corresponds to the results of B-B fitness model; q = 1.5 (dot-dashed line). 
The inset plot shows the chemical potential for another choices of q in between: q = 1.1 (open 
circles), and q = 1.3 (open triangles), respectively. Lower Panel: Fraction of the total number of 
links connected to the most connected (or "winner") node, k max /(mt), plotted as a function of T, 
shown for m = 2 and for different q-s: q = 1.0 (solid line), q = 1.1 (dashed line), q = 1.3 (dotted 
line), and q = 1.5 (dot-dashed line). 

emerge. This means the network experiencs a topological transition. In the whole condensate 
phase region (or T < Tc), the power- law solution of the degree distribution hold in FGR 
region breaks down and hence P(k) deviates from power-law distribution. However, it is 
still controlled by the nonextensive parameter q, as we can clearly see from the figure. 
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FIG. 2: Similar to Fig. [TJ but simulated for energy level distribution g(e) ~ yfi and with different 
q- value. Upper Panel: the chemical potential fi vs. temperature T with q = 1.0 (solid line), 
and q = 1.5 (dot-dashed line). The inset plot shows q = 1.1 (open circles), and q = 1.3 (open 
triangles), respectively. Lower Panel: the share of links, k max /(mt) plotted as a function of T, 
shown for m = 2 and for different q- value: q = 1.0 (solid line), q = 1.1 (dashed line), g = 1.3 
(dotted line), and q = 1.5 (dot-dashed line), respectively. 

IV. CONCLUSIONS 

In conclusion, we consider Bianconi-Barabasi fitness model of growing networks in the 
framework of non-extensive statistical mechanics via Tsallis mapping in the fitness parameter 
T]i. We find that the phenomenological nonextensive parameter q plays an important role 
in the dynamical evolution of complex networks: the underlying evolving networks may 
undergo some FGR-like phase in high temperature region with a nonextensive parameter 
controlled power-law degree distribution and it goes back to the original FGR phase of B-B 
fitness model in q — > 1 limit; on the other hand, in low temperature regime, particularly 
for T < Tc (where Tq is the critical temperature), a condensate phase emerges. For this 
condensate phase, its connectivity degree distribution P(k) deviates from the power-law 
behavior due to the fact that the power-law solution hold in FGR region breaks down, but 
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FIG. 3: The probability distribution of the degree of connectivity P(k) with q = 1.0, 1.3, and 1.5, 
for energy level distribution g{e) ~ e, at the different temperatures (from low to high). Upper-left 
Panel: at temperature T = 0.04, and q = 1.0 (open sqaures), q = 1.3 (open circles), and q = 1.5 
(open triangles); Upper-right Panel: similar to the upper-left panel but at temperature T = 0.2; 
Lower-left and -right Panels: similar to upper panels but at temperatures T = 0.8 and T = 4.0, 
respectively. 

its specific form is controlled by the nonextensive parameter q. The corresponding critical 
transition temperature T c is also identified and found decrease comparing to that of the 
original B-B fitness model. 

Tsallis formulation of nonextensive statistical mechanics not only provides a natural 
framework for describing the dynamical evolution of the growing complex networks, but 
also introduces a new complexity or entropy measure. This entropy measure is expected to 
be different from the traditional Shannon entropy or Gibbs entropy measures, e.g., those 
introduced in Refs. 32^ 36] in the context of complex networks. The current research opens 
a different angle to explore such measure, and our results may stimulate some further inves- 
tigations on this regard [371 ] . 
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